|— Anyons in small atomic Bose-Einstein condensates 
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We discuss a way of creating, manipulating and detecting anyons in rotating Bose-Einstein con- 
densates consisting of a small number of atoms . By achieving a quasidegeneracy in the atomic 
motional states we drive the system into a i-Laughlin state for fractional quantum Hall bosons. 
Localized i-quasiholes can be created by focusing lasers at the desired positions. We show how to 
manipulate these quasiholes in order to probe directly their ^-statistics. 
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The experimental achievement of Bose-Einstein con- 
densation of weakly interacting atomic gases promises 
new possibilities to study the quantum properties of 
many-body systems. As compared to other systems, 
quantum degenerate atomic gases can be easily controlled 
and manipulated by electromagnetic fields, which makes 
them ideal candidates for the study of several intrinsically 
quantum phenomena. Until now, interest has focused 
mainly on using atomic condensates to study single par- 
ticle quantum phenomena (as those occurring in atomic 
interferometry or atom optics), since for low tempera- 
tures atoms in condensates occupy essentially the same 
single particle state. In contrast, the possibility of ob- 
serving entanglement, one of the most fascinating fea- 
tures of quantum mechanics, remains almost unexplored, 
requiring a way of making the gas to effectively behave in 
a strongly interacting manner. First steps in this direc- 
tion are the recent proposals to entangle atomic beams 
and atomic spin squeezing in Bose-Einstein condensates 
(BEC) , and the numerical prediction of new correlated 
phases appearing in rotating condensates . 

In this Letter we show how a quantum degenerate Bose 
gas consisting of a small number of atoms can be used to 
study a quantum phenomenon highly collective in nature, 
namely the formation of quasiparticles exhibiting frac- 
tional statistics. This system offers the novel possibility 
of creating and manipulating anyons in a well controlled 
way that may allow for the experimental probing of their 
fractional statistics. The idea is to achieve a quasidegen- 
eracy of the atomic motional states by rotating the trap 
that confines the atoms, as it has successfully been done 
for the creation of vortices [Q. Under appropriate con- 
ditions (like two dimensional confinement), the situation 
we describe can be understood in terms of the fractional 
quantum Hall effect js) for bosons |^ . As in the theory of 
the fractional quantum Hall for electrons , elementary 
excitations exhibiting fractional statistics appear. We 
consider a situation where the atomic system is first pre- 
pared in a ^-Laughlin state, a highly correlated quantum 
liquid with nearly uniform density. We then show that 
piericing the system by offresonant laser light a single 
Laughlin quasihole localized at some chosen position zq 



can be created. This excitation involves a density pro- 
file in which exactly ^-atom has been removed at zq- In 
addition, lasers provide us with a tool for creating states 
with two quasiholes at the desired positions, and more- 
over, for creating superpositions of states having one and 
two quasiholes. Driving the superposition state along the 
proper path, as in a Ramsey-type interferometer, will al- 
low to test the fractional statistical phase. Apart from 
interest of detecting such a phase directly, the ideas de- 
veloped in the present work may pave the way for a phys- 
ical implementation for quantum information processing 
based on anyons; as reported in the most robust way 
of performing quantum computations seems to be based 
on excitations with fractional statistics since they have 
several fault-tolerant properties built-in. On the other 
hand, the conditions required to observe the effects pre- 
dicted here are very demanding, and we expect that can 
be reached in the near future only with a small number 
of atoms. 

We consider a set of N bosonic atoms confined in a po- 
tential which rotates in the y plane at a frequency Q. 
We will assume that the confinement in the z direction is 
sufficiently strong so that we can ignore the excitations 
in that direction, and we will consider that the potential 
in two dimensions is isotropic and harmonic. The Hamil- 
tonian describing this situation in a frame rotating with 
the trap is: 
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with Liz being the z component of the angular momen- 
tum of the i-th atom, and where we have used the trap 
energy, fuj, as the unit of energy, and £ — (fi/mcj)^/^ 
as the unit of length. The atoms are interacting via an 
effective contact potential, and the interacting coupling 
constant ry is related to the s-wave scattering length, a, 
and to the localization length in the z direction, £z, by 

As pointed out recently in the limit = lj the 
Hamiltonian (|l|) is formally identical to the Hamiltonian 
of electrons in the quantum Hall effect ||^ with uj playing 
the role of the magnetic field (with cyclotron frequency 
ujc — 2uj), and the usual Coulomb interaction between 
electrons replaced by a contact interaction. Regardless 
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that we are deaUng with bosons instead of fermions the 
two problems are formally identical in this limit. We will 
now study how to create and manipulate anyons in 
our twin atomic system. 

We begin by writing the Hamiltonian (|l|) in the follow- 
ing form: 

H = Hb+Hl + V. (2) 

Here, Hb = Eili -V,^/2-t-r2/2-L,2 is the quantum Hall 
single particle Hamiltonian, whose single particle energy 
levels are the Landau levels equally spaced by the cy- 
clotron energy 2hijj. The Hamiltonian Hl — (1 — Vl/uj)Lz 
is proportional to the z component of the total angular 

momentum, Lz = J^iLi^^^^ ^^'^ ^ is the interaction 
term. From now on we consider the limit in which the 
energy scales characterizing Hamiltonians Hb and V are 
much larger than the one corresponding to Hj^. This 
means that both the trap energy and the typical interac- 
tion energy are large compared to the angular momentum 
term. In this limit the ground state and elementary exci- 
tations of the system will lie on the subspace of common 
zero energy eigenstates of Hb and V . We now derive the 
spatial form of the many-body wave functions ^E" lying 
within this subspace. In order to be a zero energy eigen- 
state of Hb, ^ must lie within the subspace generated 
by tensor products of lowest Landau level single particle 
states 0, 

vI/[z]=P(zi,... ,z^)^e-l^''l'/^ (3) 

k 

where P\z\ is a polynomial in each of the atomic coordi- 
nates Zfe = a;fe + ij/fc. Let's assume that VPfz] is also an 
eigenstate of V with eigenvalue zero, and let's choose any 
pair of particles i and j. The dependence of ^'[z] on Zi 
and Zj can be reexpressed in terms of the relative and 
center mass coordinates, Zy , Zij, so that we can expand 
P\A = Em ^ijP^^ where depends on Zij and on the 
positions of all the other particles. As we are dealing 
with bosons only even values of m appear in the sum. In 
order for ^'[z] to be annihilated by the hard-core inter- 
action V , Fq must be identically zero. It follows that zfj 
is a factor of P[z] so that 

P[z]^Q[z]l[{z,-z,)\ (4) 

i<j 

We then diagonalize Hamiltonian Hl within the trun- 
cated Hilbert space of wave functions of the form speci- 
fied by (^) and (^. We note that when P[z] is a homo- 
geneous polynomial in [z] the state ^ is an eigenstate of 
Hl with eigenvalue E{M) — (1 — r2/a;)M, where the total 
angular momentum M equals the homogeneous degree of 
P[z\. It follows that the ground state of the system is the 
state with the lowest angular momentum, that is, the one 
with Q[z\ = l: 

i<j k 



E 




5 10 15 20 25 30 

M 



FIG. 1: Eigenvalues of the Hamiltonian for N — 5 bosons as a 
function of the total angular momentum M. For illustrating 
purposes we chose (1 — Q,/Lo)/r] = 0.001. Energy is measured 
in units of r]huj. 

This state is the bosonic variant of the Laughlin wave 
function for quantum Hall electrons |^ . 

We have confirmed the above arguments by diagonal- 
izing the Hamiltonian (|^) exactly for = 5 bosons. 
Fig. (0) shows the energy spectrum. There is a branch 
of states well separated in energy from the rest of the 
spectrum. These states are polynomial states of the 
form given by (||) and (P, and the ground state is the 
five atoms Laughlin state (|^), with angular momentum 
Mo = N{N -l)/2 = 20. 

To create a single fractional quasiparticles in the 
Laughlin state we insert a laser localized (within an area 
^ at some position zq. We require \zq\ to be within 
the size of the Laughlin state ( ~ 2y/N — 1). The pres- 
ence of the laser can be described by a localized repulsive 
potential, so that the new Hamiltonian of the system can 
be approximated by: 

H,^Hl + Vo J2Siz,-zo). (6) 

i 

We have studied the time evolution of the ground state 
as the intensity Vq of the laser increases with time, under 
the assumption that the system remains always in the 
truncated Hilbert space of polynomial wave functions of 
the form (^) . We note that the total angular momentum 
L does not commute anymore with the Hamiltonian (|^), 
since the (^-potential breaks the rotational symmetry. We 
have solved the dynamics exactly for the model system 
of iV = 5 bosons. Fig. || shows the time evolution of the 
ground state. For low intensities of the laser the angular 
momentum term dominates and the system remains in 
the Laughlin state. But when the laser power becomes 
sufficiently large the system evolves to a Laughlin quasi- 
hole state with wave function 

i 

where ^p[z] is the Laughlin wave function (|^). Note that 
the state -i/'zo ^ superposition of homogeneous polyno- 
mial states with angular momenta running from Mq to 
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FIG. 2: Coefficients of the ground state for the system of 
N = 5 particles as a function of time, in the Laughlin state 
ip (filled line), and in the quasihole state ipzg (dashed line). 
The laser is localized at zq = 2.5£ (the size of the droplet is 
~ 4.2^), and its intensity increases with time as Vb = O.lt 
(time in units of tj~^). The three plots at the bottom show 
the density profile of the ground state at three steps of the 
time evolution (the size of the plane is 6£ x 6£). 



Mo + N. In this quasihole state all the particles are ex- 
pelled from the position where the laser is located, so 
that the potential energy is minimized. This gain in po- 
tential energy compensates the cost of confining energy 
due to the spreading of the gas. 

It is remarkable that for intermediate laser intensi- 
ties the ground state is approximately equal to a super- 
position of the Laughlin state and the quasihole state: 
^! ^ aip + aQipzo 1 with no other states participating sig- 
nificantly in the evolution of the system. To understand 
this behavior let us consider a state quasi-orthogonal to 
this subspace, that is a state having a quasihole localized 
at some other position ^ zq. To reach this state the 
system has to pay confining energy (because of the sys- 
tem spreading around z^), but there is no gain in poten- 
tial energy since the hole has been created at the wrong 
position. It follows that any state out of the subspace 
generated by ip and ipQ will be always much higher in en- 
ergy (no matter the intensity of the laser) , and will not 
take part in the evolution of the system. 

Based on the creation of quasiholes and superposition 
states we describe now a possible experiment to reveal 
directly the statistics of anyons in a Ramsey-type inter- 
ferometer: 

l)First step: Preparation of the initial state. We fo- 
cus a laser at position zq and increase its intensity until 
a single quasihole is created. Keeping constant the in- 
tensity of this laser we then adiabatically insert another 
laser at position zi, far enough from zq (|zi — zqI ^ 
The new Hamiltonian is: Hoo — Hl + Vq ~ 
zq) + Vi J^i^i^i ~ -^i)- Following an analogous pat- 
tern to the one shown in Fig. H for a single quasi- 



hole, the system now evolves from the one-quasihole 
state to a two-quasihole state of the form: tpzg^z-^[z] — 
Yl^{zi — zo){zi — zi) ^p[z]. The crucial point is that at a 
certain point of the evolution the system reaches a super- 
position state ^ ~ {ipzQ + 'ipzczi)- This is precisely the 
superposition we need to test the statistical angle. We 
then stop the evolution of the system at this point, and 
remove instantaneously the laser located at zi. In this 
way the state of the system remains unchanged but we 
go back to Hamiltonian Hq. Note that the state 4" is no 
longer an eigenstate of this Hamiltonian. We are now in 
a position to detect the statistical phase. 

i) Second step: Statistical phase accumulation. We adi- 
abatically move the remaining laser along a path enclos- 
ing position zi, so that the time dependent Hamiltonian 
is: H{t) = Hl + VoY.^Kz^ - <t))^ with z{t) - Zi = 
\zq — zije*''*. As follows from the adiabatic theorem 
the evolved state at the end of the process will be: 

= e-'^^^+'^W^.o + e-''^^^+^^V.„,.i, (8) 

where besides the dynamical phase each state picks up 
a Berry phase. The difference between the Berry phases 
if = ifi — ip2 oi the two states reflects the extra phase 
that the quasihole at zq picks up due to the presence of 
the other quasihole at zi, thus it is the statistical phase. 

We have simulated this experiment with our 5 bosons 
model system. In order to have a more precise detection 
of the Berry phase we have considered a regime in which 
the dynamical relative phase is much smaller than the 
relative Berry phase, so that (i?i— i?2)//3 ~ 5{l—uj/il) <C 
1. We have performed the simulation for different closed 
paths and different positions of the two quasiholes. The 
relative phase emerging from the numerical calculation is 
(f = I.OSItt, so that ^ ~ V'zo.zi- Since the closed 
loop we have performed is equivalent to two consecutive 
interchanges, this result states that, as it should be, ^- 
quasiholes pick up a phase (p = -^n when we interchange 
them. Thus, they do not behave as bosons nor fermions, 
but as anyons with statistics i. 

3) Third step: Detecting the statistical phase. Suppose 
that we come back to the point at which we stopped the 
evolution of the system. This means that we instanta- 
neously restore the laser we had removed at position zi 
and continue increasing its intensity. If we had not made 
the step 2 of the experiment, the system would be in 
the state ^E* and would evolve to the two-quasihole state 
Tpzo,zi- However, after performing step 2 the state ^' has 
changed into 5*' and the system will now evolve in a dif- 
ferent way. Fig. ^ shows how the minus sign is reflected 
in the density profile of the final state. Instead of getting 
a two hole state we go back to the state with only one 
quasihole at zq. 

We discuss now the set of conditions that a system of 
N atoms must fulfill to perform an experiment as the one 
we have described above. First of all, we have made a 
two-dimension approximation, so that we need £ <^ iz- 
We have considered also a contact interaction between 
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FIG. 3: Numerical simulation of the Ramsey-like experiment 
for N — 5 bosons. Step 1: Adiabatically inserting a laser 
at position zi we drive the system from -i/j^g (plot a) into 
^ = V'zo +V'zo.zi (plot b). Step 2: We instantaneously remove 
the laser at z\ and drive the laser at zq along a closed path 
enclosing z\. The resulting state is ~ -i/i^g — V'zo.zi (plot 
c). Step 3: We restore the laser at zi and continue increasing 
its intensity. The final state (plot d) is the single quasihole 
state V'zo- contrast, if the laser at zo makes two complete 
loops around zi (step 2'). In this case the resulting state is 
again ^ ~ -tpzo + ^zo.zi (plot c'), which evolves to the state 
with two holes (plot d'). 

the atoms and this requires a < iz 0. We have cre- 
ated a state with two Laughlin quasiholes. This imphes 
a total angular momentum M = Mq+2N and thus an an- 
gular momentum energy per particle cl ^ 2(1 — fl/u!)N. 
Remember that we have projected onto the subspace of 
common zero eigenstates of Hb and V. In order for this 
projection to be valid we then need to be much smaller 
than both the trap energy and the typical interaction en- 



ergy, so that the conditions 2(1 — Q/uj)N <C 1, rj/Air are 
required For creating the i-quasiholes we had to 

focus the lasers within a distance ~ £. For a localization 
length £ ~ 1/im this implies an upper limit for the trap 
frequency of ^ IkHz. The preparation of the superpo- 
sition state used to test the statistical angle requires to 
adiabatically increase the laser intensity. In order to es- 
timate how slowly the laser needs to be inserted for a 
system of N particles we have used a two state approxi- 
mation. Confining ourselves to a basis formed by states 
ipzo "020,211 fiiid that the width of the avoided 
crossing is A A^(l - fi/a;)(|zi|/\/]V)^. This is nor- 
mally a small quantity, quickly decreasing when the hole 
is made closer to the center of the trap. The time scale 
gives us the estimation on how slowly the laser is to 
be set on/off in order to reach the necessary adiabaticity 
requirements. Finally, the most restrictive condition is 
the temperature. In order to freeze out the excitations 
we need kT/fku ^ (1 — il/oj), which together with the 
above conditions implies kT/fku ^ 1/N,ri/N |ll3| . 

In conclusion, this letter provides a way of creating 
anyons in rotating Bose Einstein condensates, and pro- 
poses an experiment in which the fractional statistics can 
be tested. While the requirements are very demanding, 
we expect that creation and (effectively) ground state 
cooling of small ensembles of bosonic atoms to be within 
experimental reach in coming years. Development of 
these experimental techniques promises the controlled 
engineering of strongly correlated, entangled states of 
atoms, with novel applications in quantum information. 
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